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We combine a pair of independent Weyl fermions to com- 
pose a Dirac fermion on the four-dimensional Euclidean lat- 
tice. The obtained Dirac operator is antihermitian and does 
not reproduce anomaly under the usual chiral transforma- 
tion. To simulate the correct chiral anomaly, we modify the 
chiral transformation. We also show that chiral gauge theo- 
ries can be constructed nonperturbatively with exact gauge 
invariance. The formulation is based on a doubler-free lat- 
tice derivative, which is a simple matrix defined as a discrete 
Fourier transform of momentum with antiperiodic bound- 
ary conditions. Long-range fermion hopping interactions are 
truncated using the Lanczos factor. 



I. INTRODUCTION 

In the continuum Euclidean path-integral, chiral 
anomaly comes from Jacobian of the fermion measure 
[1]. The point is that Tr75 of the transformed fermion 
measure gives nonzero contribution because of the infi- 
nite degrees of freedom. On the other hand, lattice field 
theory is a framework to simulate the continuum theory 
with a finite number of lattice sites. Chiral anomaly can- 
not be reproduced on the finite lattice in the same way as 
the continuum theory. We should definitely distinguish 
finite and infinite lattice formulations and concentrate on 
reproducing chiral anomaly based on the finite lattice for 
practical numerical studies. 

The species doubling problem of the lattice fermion 
[2-8] is closely related with chiral anomaly [9-12]. Ac- 
cording to the Nielsen-Ninomiya theorem, a single Weyl 
fermion cannot exist on the lattice [8]. When formulat- 
ing lattice Dirac fermion, one has almost no choice but 
to break chiral symmetry explicitly using Wilson terms 
[4]. Otherwise, one has to give up one of the other as- 
sumptions of the theorem as seen in the literature. 

Liischer's implementation of lattice chiral symmetry 
based on the Ginsparg- Wilson relation [12] was a break- 
through [13]. The most significant feature to be stressed 
is that chiral anomaly can be devised even with a finite 
lattice. In the Liischer's formulation, new chiral transfor- 
mation is introduced and chiral anomaly is obtained from 
Jacobian of the fermion measure in a different way from 
the continuum theory. The lattice index theorem holds 
for arbitrary lattice spacing [14,13] and chiral anomaly 
agrees with the continuum result in the continuum limit 
[16-19]. The lesson to be learned is that one can modify 
the axial current in order to reproduce the correct chiral 
anomaly on the finite lattice. 

In the electroweak theory, left- and right-handed 



fermions couple to gauge fields in different ways. This 
type of formulation is called chiral gauge theory. For 
consistent construction of chiral gauge theories, gauge 
symmetry needs to be maintained at the quantum level 
(see Ref. [20] , for example) . Introduction of Wilson terms 
has trouble because the mixing of left- and right-handed 
fermions complicates discussion of gauge anomaly can- 
celation [15,21]. SLAC fermion partially simplifies the 
problem because it does not use Wilson terms [3,22]. 
However, it has not been successful due to breakdown of 
locality and Lorentz invariance associated with axial cur- 
rents in the continuum limit [9] . There have been discus- 
sions that deny these defects and the non-conservation of 
the axial current [23,24]. The defects originate in the def- 
inition of axial currents with a derivative different from 
the SLAC derivative. Careful and consistent treatment 
is necessary when discussing species doubler and chiral 
anomaly with SLAC fermion. 

In this paper, we give a method to save the SLAC 
derivative. On a finite lattice, a derivative is defined 
as a discrete Fourier transform of momentum in a sim- 
ilar way to the conventional SLAC derivative. To re- 
move doubler modes, antiperiodic boundary conditions 
are chosen for the derivative in real space. The doubler 
modes at the momentum boundary are lifted and the dis- 
persion relation of the continuum theory is reproduced 
discretely. Although long-range hopping interactions ap- 
pear, the pathology of SLAC fermion is avoided and the 
correct continuum limit is guaranteed. However, long- 
range interactions are not useful in practical numerical 
calculations because sparser Dirac operator is better for 
numerical efficiency. By using the Lanczos factor tech- 
nique demonstrated in Ref. [25], we effectively truncate 
long-range hopping interactions and improve the fermion 
propagator. Since the proposed lattice derivative does 
not use Wilson terms, left- and right-handed fermions 
are completely independent. As a result, the Dirac op- 
erator constructed with the derivative has exact chiral 
symmetry. This means that the Dirac fermion does not 
reproduce chiral anomaly because the fermion measure of 
path-integral is trivially invariant under the chiral trans- 
formation on the finite lattice. To simulate the correct 
chiral anomaly, we modify chiral transformation using 
the Neuberger's solution to the Ginsparg- Wilson relation 
[7]. As a result, we obtain an anomalous Ward identity 
for the modified chiral transformation. The zero mode 
of the axial current divergence is generally non-zero and 
gives the index theorem. The absence of chiral anomaly 
under the usual chiral transformation implies the ex- 
istence of a single anomaly-free Weyl fermion. Chiral 
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gauge theories can be constructed using the single Weyl 
fermion as a building block. 

This work is a tricky extension of SLAC fermion. In the 
conventional SLAC fermion, periodic boundary condi- 
tions and an infinite lattice are assumed, where the dou- 
bler remains as a singularity at the momentum bound- 
ary. On the other hand, our formulation is based on 
antipcriodic boundary conditions and a finite lattice and 
therefore free from species doubler. In addition, the axial 
current is defined with the modified chiral transformation 
and evaluated nonperturbatively. The Nielsen-Ninomiya 
theorem does not apply to our implementation of chiral 
anomaly because the modified chiral transformation is 
used to define a symmetry. 

This paper is organized as follows. In Sec. II, the lat- 
tice derivative is defined based on the finite lattice for- 
mulation. Locality of the derivative is discussed. Long- 
range hopping interactions are truncated using the Lanc- 
zos factor. In Sec. Ill, a modified chiral transformation 
is introduced to reproduce the correct chiral anomaly. A 
way of constructing gauge-invariant chiral gauge theories 
is given in a nonperturbative way. Sec. IV is devoted to 
summary and discussions. 

II. LATTICE DERIVATIVE 

We define lattice derivative as a discrete Fourier trans- 
form of momentum 
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where n represents lattice sites and takes integer values 
between — N/2 + 1 and N/2. The lattice size N is a finite 
even number and pi corresponds to momentum. 
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Antiperiodic boundary conditions have been chosen in 
real space, \7 n+N = — V„. 1 In Eq. (1), the summation 
can be carried out easily. 
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1 If periodic boundary conditions are chosen, there appear 
degenerate zero modes on the momentum boundary. Such un- 
physical doubler modes must be removed because they cause 
errors especially when fermion mass is small. 



where s = 2im/N. The lattice derivative (1) is doubler- 
free and reproduces discretely the dispersion relation of 
the continuum theory for arbitrary lattice spacing. 
In the large N limit, Eq. (1) becomes an integral 
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where a is a lattice spacing and x — na is a space coordi- 
nate. In the continuum limit a — > 0, we obtain the first 
order derivative of the continuum theory. 
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The lattice derivative (1) is local in the continuum limit. 
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FIG. 1. The lattice derivative V(s) = V n is plotted as a 
function of s = 2nn/N for (a) N = 50 and (b) N = 1000 in a 
region \s\ < 2ir. The function V(s) has antiperiodicity of 2tt, 
V(s + 2tt) = -V(s). 

Figure 1 plots the lattice derivative V(s) = V„ with 
s = 2im/N in a range \s\ < 2n for two lattice sizes 
N = 50 and 1000. The absolute value of V(s) is large 
around \s\ = and 2w and small around \s\ — n. The 
points \s\ = and 2n are equivalent because of antiperi- 
odicity. Therefore, V(s) around \s\ = 2ir does not mean 
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severe nonlocality. The points \s\ = n give the most long- 
range hopping. As expected from Eq. (4) , locality of the 
derivative V(s) is quite good when N = 1000. On the 
other hand, decay of the derivative V(s) is slow when the 
lattice size is small. 

We are interested in constructing a theory with bet- 
ter locality for a practical purpose. In order for Eq. (1) 
to be a useful lattice derivative, long-range hopping in- 
teractions need to be truncated. However, truncation of 
hopping interactions may cause errors. We need to find 
a systematic way to reproduce spectra effectively with 
only short-range hopping interactions. For simplicity, let 
us consider a classical action for a free massless fermion 
in one-dimensional space. 

N/2 

S = a E ^m-^m-n^n- (5) 

m,n=-N/2+l 

Fourier transforms of the one-component fermions are 

N/2 

V i=-iV/2+l 
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where antiperiodicity is assumed 

1pn+N = -tpn, Tpn+N = ~^n- 

Then we have 

S = ^2(lipi,l'(l', (6) 
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where 

N/2 
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= Pi6i,i'- (?) 
with the inverse transform pi of Eq. (1) 

N ' 2 - 1 /2nln\ 
P l = - 2 E V « sin [~N~J + ^ 1 )' VAr / 2 - ( 8 ) 

To obtain this, antiperiodicity of V„, ijj n , and ^ n has 
been used. Antiperiodicity in real space gives rise to 
periodicity in momentum space, pi+n = Pi, Ci+n = (i, 
and C; + n = Ci- Some explanations would be necessary for 
the derivation of Eq. (8). Consider the matrix contained 
in Eq. (7) 

fi — y7 — i2irlm / N +i2irV n / N 



which has periodicity, C m+N>n = C m . n+N = C rn . n . In 
Fig. 2, the matrix C m ,„ is shown schematically (see the 
upper diagram). Some examples for the indices (m,n) 
are given. The matrix elements on the dotted lines are 
not contained in Eq. (7). In the upper diagram, the 
vertices (N/2, -N/2) and (-N/2, N/2) correspond to the 
points with \s\ = 2ir in Fig. 2. Using the periodicity of 
the matrix, the triangles 1 and 2 can be moved to form 
the parallelogram (see the lower diagram). As a result, 
Eq. (7) can be evaluated by calculating the summation 
for each |m — n\, which draws a line segment parallel to 
the oblique sides of the parallelogram. The last term of 
Eq. (8) is a contribution of the term with \m — n\ = 
N/2, which corresponds to the left oblique side of the 
parallelogram. 



(0,-JV/2) 




(0,JV/2) 



(0,-JV/2) 




(0,JV/2) 



FIG. 2. The matrix C m ,n is schematically shown (the up- 
per diagram). The numbers in brackets are examples of the 
indices (m,n). The matrix elements on the dotted lines are 
not contained in Eq. (7). The two triangles 1 and 2 can be 
moved to form the parallelogram shape (the lower diagram). 

In Eq. (8), we truncate long-range terms with a pa- 
rameter N c < N/2, which represents the largest distance 
of fermion hopping. 
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The truncation can be implemented to Eq. (7). 

m , =-l]T v^e- 2 ^"/^ 2 ^'"^. (10) 



The indices m and n run from -TV/ 2 + 1 to TV/ 2. The 
prime symbol means that the summation is taken over 
\m-n\ < N c , m- n + N < N c , and m-n- N > -N c . 
Fermion hopping has been restricted to a finite range. 

Figure 3 plots pi of Eq. (9) as a function of 2nl/N for 
N c = 5, 10, and 25 with a lattice size N = 50. N c = 25 
gives the exact result with no truncation, which satisfies 
the dispersion relation of the continuum theory. If one 
does not mind inclusion of long-range hopping, doubler- 
free formulation of a single Weyl fermion is possible main- 
taining the correct dispersion relation. When N c is small, 
there is also no doubler modes because of antiperiodic 
boundary conditions. Although some modes around the 
momentum boundary deviate from the correct disper- 
sion, those are not so harmful because there is no genuine 
degeneracy with low-lying modes. However, pi oscillates 
around the exact result. The oscillation tends to be large 
as N c goes to small. 
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FIG. 3. Momentum pi of Eq. (9) is plotted as a function 
oi2nl/N for N c = 5 (squares), 10 (triangles), and 25 (circles) 
with a lattice size iV = 50. iV c = 25 gives the exact result. 
The truncated results oscillate around the exact one. 



The small oscillation around the correct dispersion 
comes from the truncated terms having large n's in Eq. 
(9). As shown in Ref. [25], such oscillation can be re- 
moved by introducing the Lanczos factor, which is used 
in Fourier analysis to cancel the Gibbs phenomenon. [26] 
We modify Eq. (9) as follows: 



Pi 



n=l 



(ii) 



where 



F n = sin 



irn \ N c + 1 
is the Lanczos factor. As a result, Eq. (10) becomes 
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The final form of the action is given as 
' _ i 
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FIG. 4. Momentum pi of Eq. (11) is plotted as a function 
oi2nl/N for N c = 5 (squares) and 10 (triangles) with a lattice 
size N = 50. The result for N c = 25 (circles) is the exact one 
shown in Fig. 3. The oscillation around the exact result has 
been removed with the Lanczos factor. 
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Figure 4 plots pi of Eq. (11) improved with the Lanczos 
factor as a function of 2-kI/N for A c — 5, 10, and 25 
with TV = 50, which are compared with the exact result 
(N c = 25) shown in Fig. 3. As before, there is no doublcr 
for every A c . In addition to this, the oscillation has been 
removed with the Lanczos factor. As A c goes to large, 
the deviation around the momentum boundary tends to 
be small. In this way, we can construct a doubler-free 
ultralocal derivative. If the Lanczos factor is introduced, 
ultralocal formulation of a single Weyl fermion is possible 
maintaining almost correct dispersion relation. 



III. CHIRAL AND GAUGE ANOMALY 

On the finite four-dimensional Euclidean lattice, con- 
sider an effective action T[U] 



where 



(15) 



(16) 



is a classical action for a massless Dirac fermion coupled 
to gauge. If) = J ii Dp is a Dirac operator and the Eu- 
clidean Dirac matrices satisfy 7^ = jp and {7^,7,,} = 
28,1V. Chirality is defined with 75 = 71727374- The 
fermion variables r/i m and ip n are Grassmann valued. The 
indices m and n are four-component numbers to indicate 
lattice sites and each component runs from — A/2 + 1 to 
A/2 as before. The lattice covariant derivative 



(Dp)r 
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(17) 



v=\{v±ii) 



is diagonal with respect to the space-time indices to and 
n except for the /i-th ones. The Dirac operator D is 
antihermitian and satisfies {^,75} = 0. The classical 
action is invariant under the usual chiral transformation. 
The derivative V„ can be replaced with the truncated one 
in the same way as Eq. (14), if ultralocal construction is 
preferred. The gauge variable U m ,n(n) is a product of all 
link variables that compose a line segment between the 
two sites to and n parallel to the [i-th direction. When 
connecting the two sites with link variables along the /x- 
th direction, there arc two ways because of periodicity of 
the action. The most natural choice is the shorter path. 
One of the two ways is chosen depending on the distance 
between two sites. When \rrip — rip\ < A/2, 

T j ( \ J U m m -(i ■ ■ -Un+fr n (iTlp > Up) , . 

U m , n (n) = < jj jj , , v , (18) 

which corresponds to the hexagon sandwiched between 
the triangles 1 and 2 in Fig. 2. fi is a unit vector in the 
/U-th direction. When \m,p — Tip\ > A/2, 



U m ,m+fi ■ ■ ■ C^n+(AT-l)/i,n+Ar/i (m-p > ™p) 



U m .m — {i • • • ^n—{N+\)p^r, 



-Np (mp < Up) 



(19) 



which corresponds to the triangles 1 and 2 in Fig. 2 
and therefore intersects the boundary. The link vari- 
ables U n+ p tn are elements of a gauge group and satisfy 



U ni n+p — U\ + ^ n and U n+ p+Np, n +Nfi — U n+ p. n . In the 
continuum limit, the lattice covariant derivative becomes 



\vmADp) r 

a— »0 



iH\x-y)D^ 



where x = ma and y — na, and Dp = dp — igAp with a 
parameterization U m +p_ m = e la9A ^ x \ Eq. (17) reduces 
to the covariant derivative of the continuum theory in 
the continuum limit. 

With our Dirac operator, the usual chiral transforma- 
tion does not reproduce chiral anomaly because the clas- 
sical action and the fermion measure is invariant. To sim- 
ulate the correct chiral anomaly on the finite lattice, we 
introduce the modified chiral transformation with 9 n <^ 1 
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(21) 
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The operator G is the Neuberger's solution [7] to the 
Ginsparg- Wilson relation 75G + G75 = aG^G and has 
nothing to do with the Dirac operator (17). The mod- 
ified chiral transformation is a symmetry of the effec- 
tive action T for arbitrary lattice spacing independent of 
whether the transformation is local or global. When the 
transformation is global 9 n = 9, it is also a symmetry 
of the classical action S in the continuum limit a — > 
because the variation of the Lagrangian density induced 
by the transformation is proportional to lattice spacing. 



(24) 



Although the variation is classically zero in the contin- 
uum limit, the vacuum expectation value of the variation 
gives the index theorem for arbitrary lattice spacing. 

The axial current divergence is defined as a variation of 
the classical action under the local chiral transformation 
with 75. 

(dpJ 5 p) m = E 

n 1 ( "75)721 ,m^Pm,n2 



ni ,n 2 



(75) rn , n i ip ni ■ (25) 
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The axial current is obtained by inverting the derivative 

4 

(AO 1pn 2 



+V>n 2 (A0 
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(26) 



which is gauge invariant. In the free theory, the currents 
are local in the continuum limit because the Leibniz rule 
holds in the limit. 2 The fcrmion measure transforms as 

Dip'Dtp' = exp ^-2i QnA^ DipDtp. (27) 



Chiral anomaly 



An = tr(7 5 )„ ; . 



(28) 



is a gauge- invariant quantity. As shown in Ref. [13], the 
index theorem holds for arbitrary lattice spacing. 



y^ j A n = Indcx(G). 



(29) 



(See also Ref. [27,28].) Since the effective action L is in- 
variant under transformation of the integration variables, 



the modified axial current jf. does not conserve. 



(30) 



The Ward identity for the modified chiral transformation 
holds also for the zero mode, which gives the index the- 
orem and corresponds to the variation under the global 
transformation (24). 3 In the continuum limit, the chiral 
anomaly agrees with the continuum result [16-19] 
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where e 123 4 = +1 and F uu = [D^\ dI c) ]. The effect of 
chiral anomaly can be implemented to physical quantities 
via the modified anomalous axial current. 

Construction of anomaly-free chiral gauge theory is 
easy if our lattice derivative is used. In our formulation, 
Weyl fermions are defined with the ordinary 75 (not 75). 



2 In Ref. [9] for the conventional SLAC derivative, axial cur- 
rents are defined by introducing a derivative independent of 
the SLAC derivative, which is the cause of breaking of locality 
and Lorentz invariance . 

3 The explicit breaking term (24) is necessary for the ex- 
istence of the non-vanishing zero mode of the axial current 
divergence. 



As a result, the fermion measure of the path integral do 
not depend on gauge variables. The Weyl fermions are 
free from gauge anomaly beforehand. Consider a classical 
action for a Dirac fcrmion 



(32) 



where only the right-handed fermion is coupled to gauge. 
The left-handed fermion is redundant and does not con- 
tribute to the effective action. Once the absence of gauge 
anomaly is confirmed, the left-handed fermion can be in- 
tegrated out. The definition of the Dirac operator p is 
same as Eq. (17) except that gauge is defined as a prod- 
uct of link variables 



n-\-jj,,n 



(33) 



where A^ n = A a a n T a and P ± = (1 ± 75 )/2. p is no 
longer antihermitian. Under local gauge transformation, 
the classical action (32) is invariant. On the finite lattice, 
infinitesimal local gauge transformation 

i>' n = {l+ie a n T a P + )i> n , ^ = ^(l-<T a P_), (34) 

does not change the fermion measure 



Dip Dip = exp 



DipDip = DipDip 



and hence the effective action. A single Weyl fermion can 
exist on the lattice without violating gauge symmetry. 4 



IV. SUMMARY AND DISCUSSIONS 

We have constructed a doubler-free covariant deriva- 
tive and an anomalous Ward identity for the modified 
chiral symmetry on the lattice. The index theorem holds 
for arbitrary lattice spacing and the dependence of chiral 
anomaly on gauge fields agrees with the continuum result 
in the continuum limit. The zero mode of the anomalous 
Ward identity gives the index theorem. In this formula- 
tion, a single Weyl fermion can exist on the lattice main- 
taining gauge symmetry. Chiral gauge theories can be 
constructed nonperturbatively using the single anomaly- 
free Weyl fcrmion as a building block. 

The proposed lattice derivative is a simple matrix that 
gives the correct dispersion. We have shown that in- 
troduction of the Lanczos factor enables us to construct 
lattice derivatives with good locality. In this case, almost 
the correct dispersion is reproduced except for deviation 



The number of Weyl fermions needs to be even when there 
exists Witten's global anomaly. [29] 
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at high momentum. It depends on fermion mass how 
large the truncation parameter 7V C should be. 

The correct chiral anomaly has been obtained without 
violating locality and Lorentz invariance in the contin- 
uum limit in spite of the existence of long-range hop- 
ping interactions. This is a result of a nonperturbative 
formulation based on the modified chiral transformation 
(or equivalently the modified axial current). To repro- 
duce chiral anomaly, Wilson terms have to be introduced 
somewhere. We have introduced them only in the modi- 
fied chiral transformation to maintain complete indepen- 
dence of left- and right-handed fermions in the classical 
action. 

When calculating physical quantities dependent on chi- 
ral anomaly such as the 7/ meson mass, 75 needs to be 
replaced with 75 in the concerned vertex operators to 
include the effect of chiral anomaly. It depends on the 
construction of the operator G how precisely the effect 
of chiral anomaly is implemented with finite lattice spac- 
ing. In lattice QCD, the modification of the axial current 
does not affect physical quantities independent of chiral 
anomaly because the axial current does not couple to 
gauge directly. 
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